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We present precision calculations of the transverse-momentum spectrum, the invariant-mass 
distribution and the total cross section for Z' production at hadron colliders. We implement 
joint rcsummation at the next-to-leading logarithmic accuracy and consistently match the 
obtained result with the pure perturbative result at the first order in the strong coupling 
constant. We confront our numerical predictions with the Monte Carlo generator PYTHIA 
and with a new implementation of extra neutral gauge bosons in the Monte Carlo generator 
MC@NLO. The impact of scale dependence is also studied. 

1 Introduction 

Despite its phenomenological success, the Standard Model (SM) of particle physics is believed 
to suffer from a variety of conceptual deficiencies. In particular, it provides no fundamental 
motivation why the strong and electroweak interactions should be described by three different 
gauge groups. Grand Unified Theories allow for a unification of these groups which is broken 
to the SM at higher scale. This yields one or several extra neutral gauge Z'-bosons exhibiting 
the existence of additional U(l) symmetries^. If the Z'-bosons couple to quarks and leptons 
not too weakly and if their mass is not too large, they will be produced at the Tevatron and 
the LHC and easily detected through their leptonic decay modes. The search for these particles 
occupies therefore an important place in the present experimental program. 

When studying the transverse-momentum (px) distribution of a Z'-boson with an invariant 
mass M, it is convenient to separate the large- and small-pr regions. For the large values of px, 
the use of the fixed-order perturbation theory is fully justified, since the perturbative series is 
controlled by a small expansion parameter, a s (M 2 ), but in the small-pr region, the coefficients 
of the perturbative expansion are enhanced by powers of large logarithmic terms, ln(M 2 /pj.). 
As a consequence, the convergence of the perturbative series is spoiled as pt — > 0. Furthermore, 



when the initial partons have just enough energy to produce a Z'-boson, the mismatch between 
virtual corrections and phase-space suppressed real-gluon emission leads also to the appearance 
of large logarithmic terms [ln(l — z)/(l — #)]+, where z = M 2 /s, s being the partonic centre- 
of-mass energy. However, the convolution of the partonic cross section with the steeply falling 
parton distributions enhances the threshold contributions even if the hadronic threshold is far 
from being reached and large corrections are expected. 

Accurate calculations of pt and invariant-mass distributions must then include soft-gluon 
resummation in order to obtain reliable perturbative predictions and properly take these logarith- 
mic terms into account. We implement the joint resummation formalism^ at the next-to-leading 
logarithmic (NLL) accuracy and consistently match the obtained result with the pure perturba- 
tive result at the first order in a s §, which allows us to resum all the logarithms simultaneously. 
In addition, we compare our predictions^ with those of the next-to-leading order (NLO) Monte 
Carlo generator MC and of the leading-order (LO) Monte Carlo generator PYTHLA© 

and study the impact of the scale dependence. 



2 Theoretical framework 

2.1 Considered Z' model 

Ten-dimensional string theories with E% x Eg gauge symmetry are anomaly-free and contains 
chiral fermions as in the SM^l. After compactification, this symmetry leads to an effective Eq 
GUT group that can be broken further to^ 

E 6 50(10) x U(l) f SU(5) x U(l) x x 17(1)^. (1) 

While the Z'-bosons related to the additional U(l)^ and U(l) x symmetries can in general mix, 
we consider in this work only a TeV-scale Z x -boson and assume the to acquire its mass at 
much higher scales, as it is naturally the case in the hierarchy of symmetry breaking of Eq. (JT]). 

2.2 Joint resummation formalism at the next-to-leading logarithmic level 

In Mellin A-space, the resummed hadronic cross section for the hard scattering process 

h a h b ^ 7 ,Z,Z' ^l + r{M,p T )+X, (2) 
where a lepton pair with invariant mass M and transverse momentum p? is produced, factorizes 
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(A, b)=J2 fa/h a (N + 1) f b/hb (N + 1) H ab (N) exp{£(A, b)}. (3) 



fa,b/h ab are the A-moments of the universal distribution functions of partons a, b inside the 
hadrons h a>b , the impact-parameter b is the variable conjugate to px through a Fourier trans- 
form and the dependence on the renormalization and factorization scales fin and fip has been 
removed for brevity. The process-dependent function TL contains all the terms due to hard virtual 
corrections and collinear radiation while the process-independent Sudakov form factor Q allows 
to resum the soft-collinear radiation, embodying the all-order dependence on the logarithms. 

Once resummation has been achieved in N- and 6-space, inverse transforms have to be 
performed in order to get back to the physical space. Since the resummed exponent contains 
singularities, the integration contours of the inverse transforms must avoid hitting any of these 
poles. The 6-integration is performed by deforming the integration contour with a diversion 
into the complex 6-space^, while the inverse Mellin transform is performed following a contour 
inspired by the Minimal Prescription^! anc i the Principal Value Resummation^. 

"Note that next-to-next-to-leading order effects have recently been investigated as well^. 




Figure 1: Transverse-momentum spectra (left), mass spectra (centre) and dependence of the total cross section 
on the common factorization/renormalization scale (right) in LO QCD (dot-dashed), NLO QCD (dashed) and 
after resummation (plain). The mass spectra have been normalized to the LO QCD prediction. 



In order to keep the full information contained in the fixed-order calculation and to avoid 
possible double-counting of the logarithmically enhanced contributions, a matching procedure 
of the NLL resummed cross section to the fixed order result is performed through the formula 
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Jq(x) is the O th -order Bessel function, d 2 ^ 1 ^ ') is the fixed-order perturbative result, d 2 cr( res ) is 
the resummed cross section and dV exp ) is its truncation to 0(a s ). 



3 Numerical results 

In Fig. [1] (left), we show the pT-spectrum (integrated over the invariant mass in the range from 
900 to 1200 GeV) for Z'-bosons with a mass of 1 TeV produced at the LHC. As expected, 
the true NLO distribution (dashed) diverges as pt — > while resummation (plain) leads to a 
smooth turnover with a maximum at around 8 GeV. As can be seen from Fig. Q] (centre), the 
NLO X-factor (dashed) is increased further by the resummed contributions (plain), even if we 
are still relatively far from the production threshold. The scale uncertainty of the total cross 
section (integrated over all transverse momenta and over the invariant mass in the range from 
900 to 1200 GeV) is shown in Fig. CD (right). The LO QCD prediction (dot-dashed) does not give 
a reliable estimate of the theoretical error, since the NLO cross section (dashed) is considerably 
larger. At NLO, the factorization scale dependence is reduced as expected, but a s (fj,ji) makes 
its appearance, so that an additional renormalization scale dependence is introduced. The total 
NLO scale dependence is reduced to 9% (vertical lines) once resummation is achieved (plain). 

We confront in Fig. [2] the resummed results (plain) to the LO predictions obtained with 
PYTHIA's "power shower" (circles) and to the NLO predictions obtained with MC@NLO's 
parton shower (stars). The shaded band around the resummed predictions corresponds to the 
theoretical uncertainty induced through the simultaneous variation of the renormalization and 
factorization scales by a factor of two around the central scale M. The p^-spectra (left) are, 
for all three calculations, no longer divergent. The PYTHIA prediction rises and falls rather 
steeply around its maximum at 3 GeV, whereas the MC@NLO and resummed predictions rise 
and fall more slowly around the peak at 8 GeV. The agreement between MC@NLO and joint 
resummation is impressive, in particular for a scale choice of M/2 (upper end of the shaded 
band). The correction factors for the mass spectra (right), which have been normalized to the 
LO QCD prediction, show only a very weak mass dependence. We have multiplied the PYTHIA 




Figure 2: Mass (left) and transverse-momentum spectra (right) in PYTHIA (circles), in MC@NLO (stars), and 
after joint resummation (plain). The mass spectra have been normalized to the LO QCD prediction, and the 
renormalization and factorization scale uncertainties in the resummed predictions are indicated as shaded bands. 



mass spectrum by hand with a global .fT-factor of 1.26. Otherwise, within the statistical error 
bars, the PYTHIA K-factor would just be unity, since the normalization of the total cross section 
is not changed by the parton shower. Again, the MC@NLO X-factor agrees almost perfectly 
with the one of resummation. 



4 Conclusions 

We have improved the theoretical predictions for the production of extra neutral gauge bosons 
at hadron colliders, which are currently based on the LO Monte Carlo generator PYTHIA, by 
implementing the Z'-bosons in the MC@NLO generator and by computing their differential 
and total cross sections in the framework of joint resummation. MC@NLO and resummation 
were found to be in excellent agreement for mass spectra and pt spectra, while the PYTHIA 
predictions show significant shortcomings both in normalization and shape. The theoretical 
uncertainties from scale variations were found to be under good control. The implementation 
of our improved predictions in terms of the new MC@NLO generator or resummed K factors in 
the analysis chains of the Tevatron and LHC experiments should be straightforward and lead 
to more precise determinations or limits of the Z'-boson masses and/or couplings. 
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